Iron Fluorescent Line Emission from Young Stellar Objects in the 

Orion Nebula 
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ABSTRACT 

We present the result of a systematic search for the iron Ka fluorescent line at ~ 6.4 keV 
among 1616 X-ray sources detected by ultra-deep Chandra observations of the Orion Nebula 
Cluster and the obscured Orion Molecular Cloud 1 population as part of the Chandra Orion 
Ultra-deep Project (COUP). Seven sources are identified to have an excess emission at ~ 6.4 keV 
among 127 control sample sources with significant counts in the 6.0-9.0 keV band. These seven 
sources are young stellar objects (YSOs) characterized by intense flare-like flux variations, thermal 
spectra, and near-infrared (NIR) counterparts. The observed equivalent widths of the line cannot 
be attributed to the fluorescence by interstellar or circumstellar matter along the line of sight. 
The X-ray spectral fits and NIR colors of the 6.4 keV sources show that these sources have 
X-ray absorption of > lxlO 22 cm -2 and NIR excess emission, which is not expected when the 
fluorescence occurs at the stellar photosphere. We therefore conclude that the iron fluorescent 
line of YSOs arises from reflection off of circumstellar disks, which are irradiated by the hard 
X-ray continuum emission of magnetic reconnection flares. 

Subject headings: scattering — X-rays: stars — stars: pre-main sequence — open clusters and associa- 
tions: individual (Orion Nebula Cluster) — ISM: individual (Orion Molecular Cloud 1) 



1. INTRODUCTION 

Contrary to most X-ray emission lines that con- 
vey information of "hot" celestial objects, fluores- 
cent emission carries information of "cold" mat- 
ter in the vicinity of bright X-ray sources. For a 
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wide range of incident X-ray spectra, iron is the 
most visible element in the K fluorescent lines from 
the photoionized matter when its cosmic abun- 
dance and fluorescence yield are taken into ac- 
count (George & Fabian 1991). The Ka fluo- 
rescent line from neutral to low-ionized irons at 
~ 6.4 keV can be distinguished from Ken lines from 
highly ionized irons (~ 6.7 keV by Fe xxv and 
^6.9 keV by Fe xxvi) using X-ray charge cou- 
pled device (CCD) spectroscopy with a resolution 
of ~ 150 eV at ~ 6 keV. The line can also pene- 
trate through a large column density up to a few 
times 10 24 cm" 2 , which gives the opacity r = 1 by 
photoelectric absorption. All these features make 
iron Ka fluorescence a unique observational tool 
to investigate the environment where the contin- 
uum X-ray sources reside. 

The fluorescent line is reported in a variety 
of X-ray emitters. The two most widely known 
classes are X-ray binaries and active galactic nuclei 
(AGNs). The detections of the line also extend to 
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giant molecular clouds (e.g., Sgr B2; Koyama et al. 
1996), massive stars (e.g., 77 Car; Corcoran et al. 
1998), supernova remnants (e.g., RCW86; Vink, 
Kaastra, Bleeker 1997), and unidentified sources 
in the Galactic center (Park et al. 2004). In r\ 
Car, Corcoran et al. (1998) attributed the line to 
the fluorescence from photoionized circumstcllar 
matter along the line of sight. 

The best studied case of this emission is the 
Sun. The 6.4 keV line was first discovered dur- 
ing solar flares (Neupert et al. 1967). Following 
studies revealed that the photoionization by flare 
X-rays is more likely to produce this line than col- 
lision ionization by accelerated particles (Doschck 
et al. 1972). Parmar et al. (1984) made a system- 
atic survey of the 6.4 keV emission among ~ 600 
solar flares detected by Solar Maximum Mission. 
Their result was consistent with the calculation 
by Bai (1979), who made a model assuming the 
6.4 keV line is iron fluorescence at the solar pho- 
tosphere. 

Recently, Imanishi et al. (2001) detected the 
line emission from a class I protostar in p Ophi- 
uchi dark cloud at ~ 145 pc (de Zeeuw et al. 1999) 
with Chandra. The protostar YLW 16A produced 
a huge flare during the observation reaching a lu- 
minosity of 1.4xl0 32 ergs s _1 (0.5-9.0 keV), dur- 
ing which a 6.4 keV line emerged in its X-ray spec- 
trum with an equivalent width (EW) of < 180 eV. 
They interpreted that the emission is from a face- 
on protostellar disk irradiated by the intense X- 
rays from a magnetic reconnection event. The 
previous ASCA detection of a ~ 6.2 keV emis- 
sion line from an embedded far-infrared source in 
the R CrA cluster (Koyama et al. 1994) might be 
another example. These detections point out the 
potential importance of this X-ray line to diagnose 
deeply embedded protostars and their surround- 
ings, although the geometry of the reflection for 
the fluorescence is not certain with only a few de- 
tections. Continuum sources are presumably flare 
X-rays, but reflectors can be interstellar matter 
along the line of sight, disks, or photosphere. Sys- 
tematic studies with more samples are mandatory 
to elucidate the general picture of this faint emis- 
sion. 

This paper has two purposes: (1) to estab- 
lish the presence of the iron fluorescent line from 
young stellar objects (YSOs) with a larger num- 



ber of detections and (2) to examine the geometry 
of reflection by comparing the 6.4 keV sample to 
a control sample not showing the line. We uti- 
lize the most sensitive X-ray observation ever ob- 
tained of a star- forming region, the Chandra Orion 
Ultra-deep Project (COUP), to conduct a system- 
atic search for sources with an excess emission at 
<~ 6.4 keV. We report the detections of seven YSOs 
with the iron fluorescent line and present the in- 
terpretation that these fluorescent lines occur fol- 
lowing irradiation of disks by intense flare X-rays. 

2. OBSERVATION & DATA REDUC- 
TION 

The COUP observations were carried out in 
2003 January for 13 consecutive days with a total 
exposure time of ~ 838 ks. ACIS (Advanced CCD 
Imaging Spectrometer; Garmire et al. 2003) on- 
board Chandra (Wcisskopf et al. 2002) was used 
to take an X-ray image covering a 17'xl7' region 
in the Orion Nebula Cluster (ONC) at a distance 
of - 450 pc (Genzel & Stutzki 1989). The image 
was centered at R. A. — 5 h 35 m 17 s and Decl. = - 
5°23'40" (J2000.0) to cover the entire ONC and 
its vicinity, including the Orion Molecular Cloud 
1 cores. The data were taken with the very faint 
telemetry mode and the timed exposure CCD op- 
eration with a frame time of 3.2 s. 

Details of the data reduction are described in 
Getman et al. (2004). Briefly, after removing non- 
X-ray events and correcting for charge transfer in- 
efficiency, a source search was performed to detect 
1616 X-ray sources. They were matched with var- 
ious optical and near-infrared (NIR) catalogs for 
identification. X-ray events of each source were ac- 
cumulated from an ~ 87% encircled energy poly- 
gon at the source position. Consistent imaging, 
spectral, and timing analyses were conducted for 
all sources. In the spectral analysis, 0.5-8.0 keV 
spectra were fit by a thin-thermal plasma model 
(mekal; Mewe et al. 1985; Liedahl et al. 1995) with 
interstellar absorption (wabs), varying X-ray lu- 
minosity (£x), plasma temperature (fceT), and 
equivalent hydrogen column density (Ah) values. 
When a large residual was found, another temper- 
ature component was added to improve the fits. 
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3. SEARCH FOR FLUORESCENT LINE 
EMISSION 

3.1. Control Sample and Event Extraction 

Among the 1616 COUP sources, we first define 
a control sample to search for the 6.4 keV emission, 
which is restricted to sufficiently bright and spec- 
trally hard sources. We set two criteria of > 100 
photons and signal-to-noise ratio of > 20 in the 
6.0-9.0 keV band to obtain a sample consisting of 
127 sources. Here, we used a maximum energy of 
9.0 keV rather than 8.0 keV used by Getman et 
al. (2004) to better determine the continuum level 
beyond the iron Ka lines at 6.4-6.9 keV. 

Forty nine of these 127 sources have a "p" flag 
in the COUP source list (Getman et al. 2004), in- 
dicating that they are possibly piled-up. We de- 
veloped a method for pile-up mitigation based on 
photon extraction in an annulus around the pile- 
up core of the point spread function 1 . We accumu- 
late X-ray events from an ~ 87% encircled energy 
polygon with the central circular core removed. 
We derive an energy-dependent correction to the 
auxiliary response file for the annular region based 
on simulations using ChaRT and MARX software 
packages 2 . The correction function is a low-order 
polynomial that does not introduce artificial fea- 
tures in the corrected effective area function. Un- 
like the method by Davis (2001) which is restricted 
to mildly piled- up sources, our procedure is suc- 
cessfully applied to both mildly and severely piled- 
up sources. Events extracted from the annular re- 
gions were confirmed to meet the criteria for the 
control sample. 

3.2. Spectral Fitting 

For the 127 hard and bright COUP sources, 44 
were successfully fit with one-temperature thermal 
models and 83 were fit with two-temperature mod- 
els, following the procedure given by Getman et al. 
(2004). We made fits in the 2.0-9.0 keV band with 
-/Vjj fixed to the values tabulated by Getman et al. 
(2004). The iron metallicity value was treated as 
a free parameter, while the metallicity of other el- 
ements was left fixed to 0.3 solar. 



1 More details are given by Getman et al. (2004) and at 
http:/ /www. astro. psu.edu/users/tsujimot/arfcorr. html. 

2 See http://asc.harvard.edu/chart/ and 

http://space.mit.edu/CXC/MARX/. 



We added onto the best-fit plasma model a 
Gaussian line component at 6.40 keV with the 
intrinsic width fixed to and a free normaliza- 
tion (-/Vkq). Here, 6.40 keV is the weighted mean 
of two unresolved Ka fluorescent lines from neu- 
tral irons at 6.391 keV and 6.404 keV (Kai and 
K«2! George & Fabian 1991). The energy shift 
of the line center due to ions at different ioniza- 
tion stages makes a negligible contribution with 
< 1% for neutral to low- ionized irons (Fe i-Fe rx; 
House 1969). The intrinsic width of these lines 
(~ 0.4 cV and ~ 3.5 cV for thermal and quan- 
tum broadenings; George & Fabian 1991) and the 
width by the blending (^12 eV) are both neg- 
ligible compared to the energy resolution of the 
detector (~ 150 cV). We derive N Ka , the EW of 
the line (EWn a ), an d the null hypothesis proba- 
bility of the F test (Pp) to examine whether an 
additional model is statistically justified. With 
a careful inspection of the result, we recognized 
spectra to have a 6.4 keV feature when they meet 
(1) the 90% lower limit of Nko exceeds and (2) 
Pf is less than 1%. This suppresses the number 
of false positives to be ~ 1. We further removed 
manually a few suspicious sources that satisfied 
these criteria but did not visually show convinc- 
ing fluorescent line emission; these sources are pile- 
up sources that have distorted spectra even in re- 
duced event accumulation regions. 

This procedure produces seven positive fluores- 
cent line detections out of the 127 sources. Ta- 
ble 1 lists the name, the COUP sequence number, 
and the best-fit values both of the Gaussian and 
thermal spectral components of the seven sources, 
while the spectra of these sources are displayed in 
Figures 1. 

3.3. Light Curves 

The light curves were constructed for all the 
control sample sources using events in the 2.0- 
9.0 keV band and were segmented into pieces with 
a constant flux using a Bayesian block algorithm 
(Scargle 1998). The details on the application of 
algorithm to the COUP data are given in (Getman 
et al. 2004). 

The light curves of the seven 6.4 keV sources 
are given in Figure 2. Time-sliced spectroscopy 
of these seven sources is not feasible because of 
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Table 1 

X-ray Fitting Result of 6.4 keV Sources 



No. a 


Name 


COUP b 


Cnts° 




EW Ka ° 




AT H f 


fc B T d/ 


Lx d/ 








Scq. No. 




(10- 7 cm" 2 s- 1 ) 


(eV) 


(io- 3 ) 


(10 22 cm" 2 ) 


(keV) 


(10 30 ergs s" 1 ) 


(10 31 ergs s- 1 ) 


1 


053509.2-053058 


331 


268 


2.6 (0.1-4.6) 


126 


5.6 


1.58 


4.6 (4.3-5.6) 


3.4 (2.7-3.9) 


17 


2 


053513.6-051954 


561 


332 


3.6 (1.1-7.3) 


130 


0.8 


1.45 


2.8 (2.6-2.9) 


5.6 (5.8-6.1) 


4.4 


3 


053514.3-052232 


621 


304 


2.5 (0.5-3.9) 


111 


6.2 


5.62 


3.2 (2.9-3.4) 


4.1 (3.4-4.8) 


1.6 


4 


053514.6-052211 


647 


272 


4.1 (1.9-6.2) 


135 


5.7 


30.9 


5.7 (4.7-8.1) 


3.3 (2.9-4.1) 


5.7 


5 


053514.6-052301 


649 


107 


1.2 (0.3-2.2) 


268 


6.4 


0.93 


3.0 (2.6-3.1) 


1.4 (1.0-1.7) 


2.1 


6+ 


053519.3-052542 


1030 


461 


17 (2.5-17) 


111 


5.9 


11.5 


9.4 (6.4-14) 


11 (9.3-12) 


45 


7 


053519.6-051326 


1040 


266 


4.3 (1.5-6.6) 


122 


6.0 


1.91 


3.3 (3.0-3.5) 


6.7 (6.2-7.2) 


40 



a Pilcd-up source for which wc applied our pile-up mitigation method is marked with a dagger. 
b Scqucnce number of the COUP sources (Getman et al. 2004). 

c Counts in the 6.0-9.0 keV range. For the pile-up source, counts in the annulus region is given. 
d The 90% uncertainty is given in parentheses. 

c The fitting result of the additional 6.4 keV Gaussian to the best-fit thin-thermal plasma model. 

f The fitting result of the thin-thermal plasma model. For two-temperature spectra, the values of the higher temperature component are given. 

s The peak X-ray luminosity during the flares. The values were derived from the maximam and average count rates among segmented blocks assuming no 
spectral change. 
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Fig. 1. — Spectra and best-fit models of 6.4 keV 
sources. The upper panels show the spectra 
(pluses) and the best-fit models (solid steps) with 
a Gaussian component shown by dashed steps. 
The 6.4 and 6.7 keV lines are indicated by solid 
and broken arrows, respectively. The abscissa is 
the energy from 4.0-8.0 keV, while the ordinate is 
the spectral intensity in the unit of cnts s _1 keV" 1 . 
The lower panels show the residual to the fit, 
where the ordinate shows \ values of each bin. 



the paucity of the signal. We also checked the 
previous two short ACIS observations of the ONC 
(Feigelson et al. 2002) and found none of these 
sources have enough counts for a similar analysis. 

3.4. Stellar Counterparts of the Control 
Sample Sources 

Most of the 1616 COUP sources, thus most of 
the control sample sources, are YSOs that be- 
long to Orion Nebula population. However, a 
significant contamination comes from background 
AGNs. We estimate no more than w 6 AGNs in 
the control sample based on the following calcula- 
tion. The threshold of the control sample of 100 
counts in the 6.0-9.0 keV band is converted to 
an X-ray flux of ~ 3.3 xlO" 14 ergs s cm" 2 (2.0- 
7.0 keV) assuming a power-law spectrum of an in- 
dex typical of AGNs (-1.7) absorbed with a col- 
umn density of 10 22 cm -2 , as a rough estimate of 
the averaged value in the COUP field. Using 

/ S \~ 10 
N(> S) = 1200 - — — -j-. : r- 

v ' \2 x 10" 15 ergs s" 1 cm" 2 / 

(1) 

from Giacconi et al. (2001) where N(> S) is the 
number of background X-ray sources per square 
degree brighter than S ergs s cm" 2 , no more 
than w 6 AGNs are expected above the 100 count 
threshold. We believe four of these AGN sources 
can be identified: COUP sequence numbers 8, 723, 
1053, and 1304 (Getman et al. 2004). They exhibit 
no flare-like flux variations, show no spectral lines, 
and have best-fit UbT above 15 keV. Unphysically 
high temperature values are often introduced by 
fitting a flat power-law spectrum with a thermal 
model. We are thus confident that the remain- 
ing 123 sources are Orion population members and 
consider them to be the control sample hereafter. 

The seven sources with the 6.4 keV feature 
clearly have YSO properties: 6.7 keV line emis- 
sion arising from the Ka line of Fe xxv (Fig. 1), 
light curves showing high amplitude flux variation 
typical of YSO flares (Fig. 2), and NIR identi- 
fications except for No. 4 (Table 2). The lack 
of a NIR counterpart for No. 4 is probably due 
to an exceptionally high column density exceed- 
ing 3xl0 23 cm" 2 (Table 1). 
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Table 2 

NIR Colors of 6.4 keV Sources 



No. 




J 


H 




L' b 


L' 


TPSC C 




Counterparts 


(mag) 


(mag) 


(mag) 


Counterparts 


(mag) 




1 


05350920-0530585* 


12.5 


10.6 


9.4 








2 


387 


11.1 


9.4 


8.3 


1007 


6.7 


18 


3 
4 


428 


14.9 


11.3 


8.6 


598 


5.7 


43 


5 


446 


12.6 


11.4 


10.8 


175* 


10.5 




6 


801 


16.0 


12.4 


9.9 


42 


7.4 


41 


7 


05351965-0513264* 


14.0 


11.5 


10.2 









a J, H, and K s magnitudes arc cither from the 2MASS all-sky survey catalog (f) or McCaugh- 
rean et al. (in prep.). 

b L' magnitudes are cither from Muench et al. (2002) or Lada ct al. (2004, %). Nos. 1 and 7 
are located outside the L'-band surveys. 

c Numbers of the protostar candidate sources listed in Lada et al. (2000). 



4. DISCUSSION 
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Fig. 2. — Light curves of 6.4 keV sources in 2.0- 
9.0 keV band. The abscissa is time from the start 
of the observation, while the ordinate is the detec- 
tor count rate of each time bin. The light curves 
are binned with 2500 or 5000 s bin^ 1 according to 
their brightness. 



We have now established that seven of 123 
hard and bright COUP stellar X-ray sources show 
6.4 keV line emission. All of them exhibit power- 
ful X-ray flares and thermal plasma line emission 
typical of YSOs. Most of them also have NIR pho- 
tosphcric counterparts. We now proceed to make 
inferences about these sources concerning the na- 
ture of these stars and the origin of the fluorescent 
line emission. 

4.1. Location of the Fluorescing Material 

The observed EW values give us a constraint on 
the geometry of the continuum source and the re- 
flector that gives rise to the fluorescence. The fol- 
lowing calculation demonstrates the improbability 
that the fluorescent line arises in circumstellar or 
interstellar matter along the line of sight of the 
continuum X-rays. 

EW values of the Ka fluorescent line is related 
to the incident hot continuum and ambient cold 
material as (Liedahl 1998) 



EWi 



Ka 



Y K 



I(E Ka ) 
AO 

47T 

HE') 
E' 



Ekc 



~I(E Ka ) 
a Fe (E')dE' 



riFe(s)ds 



(2) 



Here, we suppose that the incident X-ray with a 
spectrum of 1(E) is photoelectrical^ absorbed by 
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a reflector that subtends the angle Af2 seen from 
the continuum source. The reflector, which is as- 
sumed to be optically thin to hard X-rays, has 
an iron density of n Fc (s) along the line of ray 
(s) with the photoelectric cross section of a Fe (E). 
The energy integral begins at the iron K edge 
energy \. The Ka fluorescent line with the lu- 
minosity Ikq is subsequently emitted at the en- 
ergy Eko with a fluorescence yield of Yka- As- 
suming a thermal bremsstrahlung spectrum for 
the continuum source; 1(E) oc exp (—E/k^T) and 
a Fe (E) = (E/xy 3 crFc(x), we obtain 



EW Ka 



Y K a E Ka N^Z Fe aFe(x) 



x exp 



Eko \ / k^T 



X 



/•OO 

/ exp (— x)x~ 4 dx. 

Jx/kuT 



(3) 



The integral along s is converted to an equiv- 
alent hydrogen column density (N^ — J nn(s)ds) 
using n Fc (s) = Z Fe nn(s), where nn(s) is the hy- 
drogen density and Z Fo is the elemental abundance 
of iron relative to hydrogen in the reflector. Sub- 
stituting Y Ka = 0.34 (Kortright 2001), E Ka = 
6.40 keV, X = 7.11 keV, cr Fo (x) = 2 x lO" 20 cm" 2 
(Gullikson 2001), and Z Fc = 3 x 10~ 5 (Dappen 
2000), EWko is evaluated as 



EW Ka = a 



/AO 

= a 



V 4ir 



NL 

" MeV], (4) 



10 22 cm 



where a is a constant that has a weak dependence 
on k^T. Between 3-10 keV in which our 6.4 keV 
sample sources are distributed (Table 1), a has an 
approximate value of ~ 2.5 (1.8-3.1). The value is 
slightly smaller than the one (~ 3.3) calculated for 
Seyfert galaxies with a power-law incident spec- 
trum oiI(E) cx (E/x)~ a - 7 (e.g., Krolik & Kallman 
1987), reflecting the exponential cut-off of thermal 
spectra at the hard end. 

If the fluorescence is caused by circumstellar or 
interstellar matter along the line of sight of the 
continuum X-rays, the values would be the 
same with Ah derived from the X-ray spectral fits 
(Table 1). This can account for only a few per- 
cent of the observed EW values as Ail/An < 1 
(eq. [4]). A geometry is thus required in which 
the photoionizing X-rays suffer a localized absorp- 
tion much larger than that in the line of sight. 



This is consistent with the result by Imanishi et 
al. (2001) concerning YLW 16A, where they set an 
upper limit of ~ 20 AU upon the distance from the 
protostar to the reflector, based on the lack of de- 
tectable time lag between the flare start and the 
appearance of the fluorescent emission. 

4.2. Absorption & NIR Colors of the Flu- 
orescent X-ray Sources 

Two conceivable geometries to realize localized 
reflection are: (1) the reflection by the stellar pho- 
tosphere like in the Sun; and (2) the reflection by a 
circumstellar structure like a disk. Our examina- 
tion of the X-ray absorption and the NIR colors of 
the 6.4 keV sources in comparison with the control 
sample points to a circumstellar origin. 

Figure 3 shows a scatter plot between Ah and 
k^T. The median value of Ah amongst the 123 
control sample sources is 10 219 cm~ 2 . All of the 
6.4 keV sources show larger absorption than this. 
The preference, which has a null hypothesis prob- 
ability of ~ 0.8%, is not expected if the reflec- 
tion occurs at photosphere. It would be likely, on 
the other hand, when the reflection takes place 
at disks, because YSOs with disks are younger in 
their evolution than those without disks and tend 
to have denser circumstellar matter yielding larger 
A H . 

Using J-, H-, K s -, and L'-band magnitudes 
of the NIR counterparts for the control sample 
sources, we show color-color diagrams in Figure 4. 
The NIR identifications and magnitudes follow 
Getman et al. (2004). Among 123 YSO samples, 
109 have significant J-, H-, and i^-band detec- 
tions and 66 additionally have significant L'-band 
detections to be plotted on the diagrams. The (J— 
H)/(H-K S ) and (H-K s ) / \K S -L') diagrams (Fig. 4 
a and b) complement with each other, where the 
former covers fainter sources using shorter wave- 
lengths data while the latter is more sensitive to 
cooler circumstellar matter by using longer wave- 
lengths data. At least four 6.4 keV sources (Nos. 
1, 2, 3, and 6) out of six with NIR counterparts 
lie outside the region of reddened dwarfs and gi- 
ants indicating the presence of NIR-emitting inner 
disks. The preference to have NIR excess emission 
is consistent the idea that fluorescent reflection oc- 
curs at disks. 
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Fig. 3. — Scatter plot of plasma temperature and 
column density of the control sample (squares). 
Sources with the 6.4 keV feature are marked 
filled with their source numbers. The dashed-and- 
dotted lines indicate the median values of Nu and 
k^T of the control sample. The dotted curve is 
the column density to have the opacity r = 1 to 
the X-rays at a given temperature, which roughly 
indicates the observational bias. 
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Fig. 4.— {J-H)/{H-K S ) and {H-K S )/(K S -L') di- 
agrams (a and b) of the control sample sources. 
Squares are for sources with larger k^T and A^h 
than the medium values, while diamonds are for 
remainders. The 6.4 keV sources are marked filled 
with a label of their source numbers. The red- 
dening lines are shown with dashed lines from the 
intrinsic colors of dwarfs and giants {thick solid 
curves; Tokunaga 2000) and the classical T Tauri 
star locus {thick solid lines] Meyer, Calver, & Hil- 
lenbrand 1997). No conversion of the color sys- 
tems is made to the data. 
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4.3. Fluorescence and Flares 
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All the 123 control sample sources have more 
Bayesian block segments than two, indicating flux 
variation. Most of them exhibit variation typical 
to flares with a fast rise and slow decay profile. All 
of the seven sources with iron fluorescence have 
large amplitude flares (Fig. 2). 

Our sample is heavily biased for flare sources. 
Therefore, we do not conclude that a flare is a 
required factor for fluorescence, despite the facts 
that all fluorescence sources exhibit flares and that 
we did not favor for flare sources in identifying the 
fluorescence line by fitting time-integrated spec- 
tra. 

The flare amplitude, however, appears to be re- 
lated to the fluorescence. Figure 5 shows the his- 
tograms of flare amplitude of the control sample 
and 6.4 keV sources, where we define the flare am- 
plitude as the ratio between the maximum and 
the minimum count rates of all Bayesian block 
segments. Six out of seven 6.4 keV sources have 
larger amplitudes than the median value (^10 142 ) 
among the control sample, making the null hy- 
pothesis probability of this preference to be ~6%. 
Large amplitude flares may have a preferable ge- 
ometry for fluorescence to occur, which should be 
examined by future studies. 

4.4. Constraints on the Geometry of the 
Fluorescent Sources 

We consider that the lack of 6.4 keV detections 
from the other 116 sources is not completely at- 
tributable to a statistical inadequacy of the data. 
About 14% of the 116 YSO samples without a 
6.4 keV detection show: (1) iVn and k^T larger 
than the median values and (2) significant NIR ex- 
cess emission (squares located right to the middle 
reddening lines of Fig. 4 a and b). Most of them 
have more counts in the 6-9 keV band than the 
weakest of the seven 6.4 keV sources (107 counts 
for No. 5; Table 1), which is close to the threshold 
for the control sample. The absence of the flu- 
orescence line thus appears to be a physical not 
data-dependent statistical characteristic. 

One idea is the disk orientation, if the fluores- 
cence occurs at disks illuminated from atop (the 
lamp post geometry) . The geometry is considered 
a likely explanation for the iron fluorescent emis- 
sion from the accretion disks in Seyfert galaxies 
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Fig. 5. — Histograms of flare amplutude for the 6.4 
keV sources (solid) and the control sample sources 
(dashed). The frequency of the control sample 
sources is reduced by a factor of 10 to facilitate 
comparison. 

and X-ray binaries (Nayakshin & Kallman 2001). 
Face-on inclination is more favorable than edge- 
on inclination for the fluorescence X-rays to es- 
cape. George & Fabian (1991) presented a relation 
between the expected EW of the fluorescent line 
and the disk inclination with a radiative transfer 
treatment. The calculation assumed a configura- 
tion where a continuum source is at a height of 
O.Olr above the center of a disk with a radius of r. 
The EW value decreases drastically for large in- 
clination angles relative to the normal of the disk 
(see George & Fabian 1991, Fig. 14). It may be 
the case that the bright sources without a 6.4 keV 
detection have unfavorable inclination angles seen 
from us. 

A simple consideration on the observed EW 
values also supports this picture. The iron fluo- 
rescence can come out from any depth from the 
surface of disks corresponding to different 
values, but we speculate that the observed flu- 
orescence is dominated by the light reflected at 
iVy ~ 10 24 cm~ 2 , where 6.4 keV X-rays suffer pho- 
toelectric absorption of the opacity r < 1. We 
would expect a smaller EW-p e if 1S smaller 
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than ~ 10 24 cm~ 2 by equation (4). The claim 
is also the case if is larger, because the emit- 
ted 6.4 keV X-rays would be severely attenuated 
through the reflector (Makishima 1986). By sub- 
stituting N{i ~ 10 24 cm" 2 and EW Fe ~ 140 eV, 
which is the mean value of all seven sources (Ta- 
ble 1), we can derive (Af2/47r) ~ 0.5. This is 
consistent with the assumption that the disk is il- 
luminated in a lamp post configuration to realize 
the subtended angle, and the disk is oriented face- 
on so that the fluorescent X-rays suffer no extra 
attenuation through the disk. 

The detection of the 6.4 keV line has an astro- 
physical implication on X-ray ionization of disks. 
Many studies (e.g., Glassgold, Feigelson, & Mont- 
merle 1999) discussed that X-ray photoionization, 
in addition to cosmic-ray collision ionization, can 
be a dominant ionization source for circumstellar 
material around forming stars, which is an impor- 
tant process to couple the gas with the ambient 
magnetic field. Our detection of the recombina- 
tion line at 6.4 keV is a direct observational proof 
that the X-ray photoionization does take place at 
disks. 

Finally, we note that the reflection may occur 
in different disk geometries such as X-ray irradi- 
ation at the inner edge of the disk, or at other 
circumstellar structures such as funnel flow, jet, 
or wind columns; these structures are associated 
with sources at an early stage of YSO evolution. 
So far, we have no observational argument to dis- 
tinguish among these possibilities. The measure- 
ment of the delay time from the onset of flares to 
the appearance of the 6.4 keV line, as was tried 
by Imanishi et al. (2001) for an upper limit value, 
will bring a key observational basis by future X- 
ray experiments like Constellation- X. 

The authors express gratitude for Konstantin 
Getman for his help in data reduction, and Sal- 
vatore Sciortino and Koji Mori for improving 
the manuscript. COUP is supported by Chan- 
dra guest observer grant SAO GO3-4009A (E. 
D. Feigelson, PI) and the ACIS Team contract 
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the Promotion of Science. 
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